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ABSTRACT 



We analyzed the available LIGO data coincident with GRB 070201, a short duration hard spectrum 7-ray 
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burst whose electromagnetically determined sky position is coincident with the spiral arms of the Andromeda 
galaxy (M31). Possible progenitors of such short hard GRBs include mergers of neutron stars or a neutron star 
and black hole, or soft 7-ray repeater (SGR) flares. These events can be accompanied by gravitational-wave 
emission. No plausible gravitational wave candidates were found within a 180 s long window around the time 
of ORB 070201. This result implies that a compact binary progenitor of GRB 070201, with masses in the range 
1 Mq < mi < 3 Mq and 1 Mq < m2 < 40 M©, located in M31 is excluded at > 99% confidence. Indeed, if 
GRB 070201 were caused by a binary neutron star merger, we find that Z) < 3.5 Mpc is excluded, assuming 
random inclination, at 90% confidence. The result also implies that an unmodeled gravitational wave burst 
from GRB 070201 most probably emitted less than 4.4 x 1O~'^M0C^ (7.9 x 10^^ ergs) in any 100 ms long period 
within the signal region if the source was in M31 and radiated isotropically at the same frequency as LIGO's 
peak sensitivity {f ^ 150 Hz). This upper limit does not exclude current models of SGRs at the M31 distance. 

Subject headings: gamma-ray bursts (GRB 070201) - gravitational waves - compact object mergers - soft 
gamma-ray repeaters 
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1. INTRODUCTION 

Gamma-ray bursts (GRBs) are intense flashes of 7-ray s 
which are observed to be isotropicaUy distributed over the 
sky (see, e.g.: Klebesadel et al. 1973; Piran 2005; Meszaros 
2002, and references therein). The variability of the bursts on 
short time scales indicates that the sources are very compact. 
Combined observations, using 7-ray and x-ray satellites such 
as Vela, CGRO, BeppoSax, HETE, Swift, Konus-Wind, and 
INTEGRAL (see Klebesadel et al. 1973; Meegan et al. 1992; 
Paciesas et al. 1999; Frontera et al. 2000; Mazets & Golenet- 
skii 1981; Gehrels et al. 2004, and references therein), mis- 
sions, as well as by the Interplanetary Network (IPN), with 
follow-up by X-ray, optical and radio telescopes of the region 
around GRBs, have yielded direct observations of afterglows 
from ?^350 GRBs. In turn, host galaxies were identified for 
many GRBs and redshifts were determined for ^\25 bursts. 
The redshifts indicated that GRBs are of extra-galactic origin. 
Two types of GRBs are distinguished by their characteristic 
duration (see Kouveliotou et al. 1993; Gehrels et al. 2006) 
and are understood to have different origins. 

Long GRBs have duration > 2 s. Detailed observations of 
long GRBs demonstrate their association with star-forming 
galaxies ranging up to a redshift of z 6.3 (see Kawai et al. 
2006; Watson et al. 2006; Jakobsson et al. 2006, and ref- 
erences therein). Furthermore, several nearby long GRBs 
have been spatially and temporally coincident with super- 
novae (e.g. Campana et al. 2006; Malesani et al. 2004; Hjorth 
et al. 2003; Galama et al. 1998; Woosley & Bloom 2006). 

Short GRBs have duration ^ 2 s. The progenitors of short 
GRBs are not so well understood. While there are associa- 
tions with distant galaxies of different types and different star 
formation histories, there are also powerful bursts of 7-ray s 
from Galactic sources, such as SGR 1806-20 (Nakar et al. 
2006; Hurley et al. 2005; Palmer et al. 2005). However, statis- 
tical analyses indicate that at most 15% of known short GRBs 
can be accounted for as soft 7-ray repeaters (SGRs) (Nakar 
et al. 2006; Chapman et al. 2007). Moreover, the spectral 
characteristics and energetics of some observed short GRBs 
and their afterglows seem to contradict this hypothesis in most 
cases (Nakar et al. 2006). The current leading hypothesis to 
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explain most short GRBs is the merger of neutron star or neu- 
tron star + black hole binaries (see for example Nakar 2007; 
Bloom et al. 2007, and references therein). However, to date 
no observations have definitively confirmed the association 
between short GRBs and binary mergers. 

Therefore, given the candidate sources, it is plausible that 
GRB central engines are also strong gravitational wave (GW) 
emitters at frequencies accessible to ground-based detectors 
like LIGO, GEO-600, and Virgo (Abbott et al. 2005a; Acer- 
nese et al. 2006; Willke et al. 2002; Kochanek & Piran 1993; 
Abbott et al. 2006a, 2005b; Finn et al. 2004). Bursts of grav- 
itational waves are expected to be emitted during the GRB 
event, with a characteristic duration comparable to that of the 
associated GRB, though the amplitude and frequency spec- 
trum of the gravitational- wave burst are unknown. In the case 
of short GRBs produced by compact binary mergers, grav- 
itational waves with relatively well-modeled amplitude and 
frequency evolution will be emitted during the inspiral phase 
of the binary system, preceding the event that produces the 
GRB. 

GRB 070201 was an intense, short duration, hard spec- 
trum GRB, which was detected and localized by 4 IPN 
spacecraft (Konus-Wind, INTEGRAL and MESSENGER); it 
was also observed by Swift (BAT) but with a high-intensity 
background as the satellite was entering the South Atlantic 
Anomaly (Golenetskii et al. 2007b). The burst light-curve 
exhibited a multi-peaked pulse with duration ^ 0.15 s, fol- 
lowed by a much weaker, softer pulse that lasted ^ 0.08 s. 
Using early reports, Perley and Bloom (Perley & Bloom 
2007) pointed out that the initial IPN location annulus of 
the event intersected the outer spiral arms of the Andromeda 
galaxy (M31). The refined error box, centered ^ 1.1 de- 
grees from the center of M31, was later reported (Pal' Shin 
2007; Hurley et al. 2007), and it still overlaps with the spi- 
ral arms of M31 [see Figure 1 and (Mazets et al. 2007; IPN3 
2007)]. Based on the Konus-Wind observations (Mazets et al. 
2007; Golenetskii et al. 2007a), the burst had a fluence of 
1.57(-0.21,-h0.06) X 10-5 erg.cm-2 in the 20 keV - 1 MeV 
range. 

It was also pointed out (Golenetskii et al. 2007a) that if the 
burst source were actually located in M31 (at a distance of 

770 kpc) the isotropic energy release would be ~ 10"^^ erg, 
comparable to the energy release in giant flares of soft 7-ray 
repeaters: e.g., the 5* March 1979 event from SGR 0526-66 
(~ 2 X 10"^"^ erg in the initial pulse) and the 27* December 
2004 event from SGR 1806-20 (- 2 x 10^^ erg). Conversely 
if the event had an isotropic energy release more typical of 
short hard GRBs, e.g., - 10"^^ - 10^^ erg (Berger 2007), then 
it would have to be located at least ~ 30 times further than 
M31 (i.e., further than 23 Mpc). 

At the time of GRB 070201, the Hanford detectors 
of the Laser Interferometer Gravitational-wave Observatory 
(LIGO) (Abbott et al. 2004) were stable and recording 
science-quality data, while the LIGO Livingston, GEO-600, 
and Virgo detectors were not taking data. The LIGO data 
around GRB 070201 were searched for evidence of a gravita- 
tional wave signal from compact binary inspiral or the central 
engine of the GRB itself. 

A standard measure of the sensitivity of a detector to grav- 
itational waves is the distance to which an optimally oriented 
and located double neutron star binary would produce a re- 
sponse in the datastream that, when optimally filtered for the 
inspiral waves, peaks at a signal to noise ratio of 8 (see, e.g. 
Abbott et al. 2005b, and references therein). At the time of 



Fig. 1 .— The IPN3 (IPN3 2007) (7-ray) error box overlaps with the spiral 
arms of the Andromeda galaxy (M31). The inset image shows the full error 
box superimposed on an SDSS (Adelman-McCarthy et al. 2006; SDSS 2007) 
image of M3 1 . The main figure shows the overlap of the error box and the 
spiral arms of M31 in UV light (Thilker et al. 2005). 

GRB 070201, this distance was 35.7 Mpc and 15.3 Mpc for 
the Hanford 4 km and 2 km detectors, respectively. However, 
the sensitivity of a detector to a gravitational wave depends on 
the location of the source on the sky and on the polarization 
angle of the waves. In the case of compact binaries, it also 
depends on the inclination angle of the orbital plane relative 
to the line of sight. At the time of GRB 070201, the binary 
inspiral reach in the direction of M31 was only about 43% of 
this maximum. More details of the instrumental sensitivity 
can be found in Sec. 2. 

The search for gravitational waves from a compact bi- 
nary inspiral focused on objects with masses in the ranges 
1 Mq < mi < 3 Mq and 1 Mq < m2 < 40 M©. The core 
of the search is matched filtering, cross-correlating the data 
with the expected gravitational waveform for binary inspiral 
and uses methods reported previously (see, e.g Abbott et al. 
2005b, and references therein). Uncertainties in the expected 
waveforms can lead to decreased sensitivity of the search to 
the gravitational wave signal from the inspiral phase; this is 
particularly true of systems with higher masses and systems 
with substantial spin (Grandclement et al. 2003). This is ac- 
counted for by studying the dependence of sensitivity of the 
search to a variety of model waveforms based on different ap- 
proximation methods. 

The search for more generic transient gravitational waves 
coincident with the 7-ray burst is based on cross-correlating 
data from two detectors and does not make use of a specific 
model for the gravitational wave signal. This is an appro- 
priate method when the gravitational wave signal is not well 
modeled theoretically, such as signals from the actual merger 
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phase of a compact binary system or the core collapse phase 
of a supernova event. 

The remainder of the paper is organized as follows. In 
Sec. 2, we discuss the LIGO detectors and the data taken 
around the time of GRB 070201. In Sec. 3, we report on 
the inspiral gravitational wave search, briefly reviewing the 
methods and algorithms used, and concluding with the astro- 
physical implications of the search for the GRB 070201 event. 
In Sec. 4, we report on the search for burst-like gravitational 
wave signals and present the astrophy sical implications of that 
search. The software used in this analysis is available in the 
LIGO Scientific Collaboration's data analysis code archives 
(LIGO Scientific Coflaboration 2007b) Since no plausible 
gravitational wave signal was detected above the background 
either in the inspiral or the burst search, we present the astro- 
physical implications of these results on the understanding of 
short GRBs in Sec. 5. 

2. LIGO OBSERVATIONS 

LIGO is comprised of three instruments at two geographi- 
cally distinct locations (a 4 km detector and a 2 km detector 
at Hanford Observatory, referred to as HI and H2, and a 4 km 
detector at Livingston Observatory, referred to as LI). Five 
science runs have been carried out to date. GRB 070201 oc- 
curred during the most recent science run, called S5, which 
started on November 4th, 2005 and ended on October 1st, 
2007. All three LIGO detectors were operating at their design 
sensitivity (LIGO Scientific Collaboration 2007a) throughout 
the S5 run. 

The LIGO detectors use suspended mirrors at the ends of 
kilometer-scale, orthogonal arms to form a power-recycled 
Michelson interferometer with Fabry-Perot cavities. A gravi- 
tational wave induces a time-dependent strain h(t) on the de- 
tector. While acquiring scientific data, feedback to the mirror 
positions and to the laser frequency keeps the optical cavi- 
ties near resonance, so that interference in the light from the 
two arms recombining at the beam splitter depends on the dif- 
ference between the lengths of the two arms. A photodiode 
senses the light, and a digitized signal is recorded at a sam- 
pling rate of 16384 Hz. The data are calibrated and converted 
into a strain time series. 

The LIGO detectors have a sensitive frequency band ex- 
tending from ~40 Hz to ~2000 Hz, with the maximum sensi- 
tivity at ~ 150 Hz, which is limited at low frequencies by seis- 
mic noise and at high frequencies by laser shot noise. In addi- 
tion, environmental disturbances, control systems noise, and 
other well understood noise sources result in a non- stationary 
and non-Gaussian background. 

2.1. LIGO observations coincident with GRB 070201 

At the time of the GRB trigger, both LIGO Hanford de- 
tectors were stable and recording science quality data. These 
detectors had been in science mode for more than 14 hours 
before the GRB trigger, and stayed in science mode for more 
than 8 hours after the GRB trigger, providing ample data for 
background studies. 

An asymmetric 180 s on- source segment, -120/+60 s about 
the GRB trigger time, was searched for gravitational-wave 
signals. This choice (Abbott et al. 2005a, 2007a) is conserva- 
tive enough to accommodate inspiral type signals, trigger time 

The search for inspiraUng binaries (Sec 3.) used LAL and LALAPPS 
with tag s5_grb070201_20070731 and the burst search (Sec. 4) used the 
MATAPPS package grbxcorr with tag grbxcorr_rl 



ambiguities, and theoretical uncertainties. We also implicitly 
assume that the propagation speed of the gravitational waves 
is the speed of light. The significance of candidate events was 
evaluated using studies covering several hours of off-source 
data from the same science mode stretch outside of, but near 
to, the on- source segment. 

The ideal response of a detector to an incident gravitational 
wave is a weighted combination of the two underlying gravi- 
tational wave polarizations denoted by h+(t) and hx(t)'. 

h(t) = F^(0,^,^)h^(t) + Fx(0,^,^)hAt). (1) 

The dimensionless weighting amplitudes, or antenna factors, 
F+ and Fx depend on the position (6>, (/>) of the source relative 
to the detector and the gravitational wave polarization angle 
i/j. For the location of GRB 070201, the root-mean- square 
(RMS) antenna factor, Frms, for both colocated and coaligned 
Hanford detectors was 

Frms = ^{Fi + Fl)/2 = 0.43/ V2 = 0.304 , (2) 

a combination which does not depend on the polarization an- 
gle i/j. Despite the sub-optimal location of GRB 070201 for 
the LIGO Hanford detectors, they still had significant sensi- 
tivity for the polarization state compatible with the detector. 

2.2. Data quality for the times surrounding the GRB 070201 

trigger 

A suite of data quality tests were applied to LIGO data. No 
anomalous behavior was found in either instrument at the time 
of GRB 070201. On the other hand, a number of data quality 
issues were identified in the off-source time used for back- 
ground estimation (which amounted to 60084 s, or 16.7 hr). 
Triggers were excluded from 530 s of coincident, off- source 
data so identified, or 0.9% of the off- source time. 

Overflows in digital signals used in the feedback control 
systems were responsible for 29 s in HI and 29 s in H2 of 
excluded time. Seismic noise in the 3-10 Hz band known to 
produce false alarms in HI was used to veto 160 s of data. 
Disturbances that produced a loss in power in the H2 detector 
arm cavities larger than 4% were also vetoed, amounting to 
163 s, which include lis when there were overflows in H2. 
No such fluctuations in arm power were observed in HI. 

Additionally, in the search for a compact binary progenitor, 
there were losses in off-source live-time due to quantization 
on 180 s intervals. Each of these intervals was intended to be 
a trial treated the same as the central, on-source interval, for 
use in background determination. 

For the burst analysis, three hours of data were used for 
the purpose of background estimation. The same data qual- 
ity flags were considered as were used in the inspiral search 
but, due to the shorter length of the background used, only one 
data quality flag (an overflow in the H2 signal) was applied ve- 
toing one of the 180 s segments in the three-hour background 
period. 

Finally, 160 s of the off- source time were excluded from 
this data analysis, as it contained simulated signals. These 
were injected intentionally into the hardware at predetermined 
times to validate the detector response and signal detection 
algorithms. 

3. SEARCH FOR GRAVITATIONAL WAVES FROM A COMPACT 
BINARY PROGENITOR 

A number of searches for gravitational waves from com- 
pact binaries have been completed on the LIGO data (Abbott 
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et al. 2005b, 2006b, 2005c, 2007b). Similar search methods 
were applied to the on-source time around GRB 070201 (Ab- 
bott et al. 2007b). In this section we briefly describe those 
methods, report the results of the search, and discuss their in- 
terpretation. 

3.1. Search Method 

The core of the inspiral search involves correlating the 
LIGO data against the theoretical waveforms expected from 
compact binary coalescence, i.e., matched filtering the 
data (Wainstein & Zubakov 1962). The gravitational waves 
from the inspiral phase, when the binary orbit tightens under 
gravitational-wave emission prior to merger, are accurately 
modeled in the band of LIGO sensitivity for a wide range of 
binary masses (Blanchet 2006). The expected gravitational- 
wave signal, as measured by LIGO, depends on the masses 
and spins of the binary elements, as well as the spatial loca- 
tion, inclination and orientation of the orbital axis. In general, 
the power of matched filtering depends most sensitively on ac- 
curately tracking the phase evolution of the signal. The phas- 
ing of compact binary inspiral signals depends on the masses 
and spins, the time of merger, and an overall phase. In a search 
for gravitational waves from compact binaries, one therefore 
uses a discrete set of template waveforms against which the 
data are correlated. 

In this search, we adopt template waveforms which span 
a two-dimensional parameter space (one for each component 
mass) such that the maximum loss in signal-to-noise (SNR) 
for a binary with negligible spins would be 3%. While the 
spin is ignored in the template waveforms, we show below 
that the search is still sensitive to binaries with most physi- 
cally reasonable spin orientations and magnitudes with only 
moderate loss in sensitivity. To generate a GRB, at least one 
of the objects in a compact binary must be a material object, 
probably a neutron star, while the second object must either 
be a neutron star or a stellar mass black hole with low enough 
mass (Vallisneri 2000; Rantsiou et al. 2007) to cause disrup- 
tion of the neutron star before it is swallowed by the hole. 
The mass-parameter space covered by the templates is there- 
fore 1 Mq < mi < 3 Mq and 1 M© < m2 < 40 Mq. The 
number of template waveforms required to achieve this cov- 
erage depends on the detector noise curve; at the time of the 
GRB, 7171 and 5417 templates were required in HI and H2, 
respectively. 

The data from each of the LIGO instruments are filtered 
through the bank of templates. If the matched filter signal- 
to-noise exceeds a threshold p*, the template masses and the 
time of the maximum signal-to-noise are recorded. For a 
given template, threshold crossings are clustered using a slid- 
ing window equal to the duration of the template as explained 
in (Allen et al. 2005). For each trigger identified in this way, 
the coalescence phase and the effective distance — the distance 
at which an optimally oriented and located binary would give 
the observed signal-to-noise assuming masses to be those of 
the template — are also computed. Triggers identified in each 
instrument are further required to be coincident in the time 
and mass parameters between the two operating instruments 
taking into account the correlations between those parameters. 
This significantly reduces the number of background triggers 
that arise from matched filtering in each instrument indepen- 
dently. Because HI was more sensitive than H2, two different 
thresholds were used in the matched filtering step: p* =5.5 
in HI and p* = 4.0 in H2. This choice takes advantage of the 
better sensitivity in HI while still using H2 to reduce the rate 



of accidentals. 

To further reduce the background, two signal-based tests 
are applied to the data. First, a statistic (Allen 2005), which 
measures the quality of the match between the data and the 
template, is computed; triggers with large are discarded. 
Second, the veto (Rodriguez 200'^ ), which looks at the time 
the x^ statistic stays above a threshold, is applied. 

The SNR and x^ from a single detector are combined into 
an effective SNR (Abbott et al. 2007b). The effective SNRs 
from the two detectors are then added in quadrature to form a 
single quantity p^^^ which provides good separation between 
signal candidate events and background. The final list of co- 
incident triggers are then called candidate events. 

3.2. Background and Results 

Gravitational-wave detectors are susceptible to many 
sources of environmental and intrinsic noise. These sources 
often result in non-Gaussian and non- stationary noise back- 
grounds. In the case of HI and H2, which share the same 
vacuum enclosure, these backgrounds are correlated. To esti- 
mate the background in this search, an equal number of 1 80 s 
off- source segments were selected to the past and future of the 
GRB trigger. All of the data, including the on-source segment, 
were analyzed using the methods described above. Triggers 
arising from the on-source segment were then removed, as 
were triggers within bad quality segments, leading to an esti- 
mate of the number of accidental triggers per 180 s segment. 
A total off-source time of 56340 s was analyzed, correspond- 
ing to 313 trials of 180 s. The mean rate of coincidences was 
2.4 per 180 s segment. 
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Fig. 2. — A cumulative histogram of the expected number of background 
triggers in 180 s based on the analysis of the off-source times (pluses) as 
a function of the effective signal-to-noise ratio (Abbott et al. 2007b). The 
shaded region indicates the Icr variation in the background estimate observed 
in the off-source times. 



Figure 2 shows the expected number of coincidences above 
each Pgff value in 180 s based on the analysis of the off-source 
times (Abbott et al. 2007b). No candidates were observed 
in the on-source time. Therefore, no plausible gravitational 
wave signals from compact binary coalescence were identi- 
fied around the time of GRB 070201. 

3.3. Astrophysical Interpretation 

The observations reported here can be used to constrain the 
distance to the GRB assuming it was caused by a compact bi- 
nary merger. With similar considerations, one can also eval- 
uate the probability that a compact binary progenitor at the 
distance of M31 was responsible for GRB 070201. 
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We discover these bounds by computing the HkeHhood of 
our observation, namely the probabiHty that no signal would 
be observed in the on-source time, given the presence of a 
compact binary progenitor with various parameters. Denote 
the gravitational- wave signal by h{t\m2^D^ jl) where m2 is the 
mass of the companion, D is the physical distance to the bi- 
nary, and jl= {mi,ii,i2,^,<l>o,^o} is the mass of the neutron 
star, the spins, the inclination, the coalescence phase, and the 
coalescence time. The probability of interest is then 

p[0\hit-m2,D)] = j p{pi)p[0\hit-m2,D,pi)]dpi (3) 

where the nuisance parameters jl are integrated over some 
prior distribution p{jl). This integration was performed by 
injecting simulated signals into the data streams of both de- 
tectors according to the desired prior distribution, and evaluat- 
ing the efficiency for recovering those injections as candidate 
events (as described in Sec. 3.1), as a function of m2 and D. 
We choose uniform priors over mi (IM© <mi< 3Mq), ^q, to 
and the polarization angle; the priors for spin and inclination 
L are discussed below. 

Astrophysical black holes are expected to have substantial 
spin. The maximum allowed by accretion spin-up of the hole 
is (a/M) = (cS/GM^) < 0.9982 (Thorne 1974) in units of the 
Kerr spin parameter (S is the spin angular momentum of the 
black hole). More detailed simulations and recent observa- 
tions provide a broad range of values (O'Shaughnessy et al. 
2005) with a maximum observed spin (a/M) > 0.98 (Mc- 
Clintock et al. 2006). The maximum spin that a neutron 
star can have is estimated from a combination of simulations 
and observations of pulsar periods. Numerical simulations of 
rapidly spinning neutron stars give (a/M) < 0.75 (Cook et al. 
1994); the maximal spin of the observed pulsar sample may 
be substantially lower than that. In our spinning simulations, 
we adopted a distribution in which the spin magnitudes are 
uniformly distributed between zero and (a/M) = (cS/GM^) = 
0.98 and (a/M) = (cS/GM^) = 0.75 for the black holes and 
neutron stars respectively, while the direction of each spin is 
uniform over the sphere. There is strong evidence that short 
GRBs are beamed (see, e.g., Soderberg et al. 2006; Nakar 
2007; Burrows et al. 2006, and references therein), although 
probably less beamed than long bursts (Grupe et al. 2006). If 
this is the case, the most likely direction for beaming is along 
the total angular momentum vector of the system. For binaries 
with small component spins, this will correspond to the direc- 
tion orthogonal to the plane of the orbit. Hence the inclination 
angle of the binary, relative to the line of sight, is most likely 
to be close to zero. However, since zero inclination is the best 
case for detection of gravitational waves, a uniform prior on 
cos L provides a conservative constraint. We drew cos l from a 
uniform prior. 

Figure 3 shows the contours of constant probability 1 - 
p[0\ h(t',m2^D)]. Compact binaries corresponding to parame- 
ters (m2,D) in the darkest- shaded region are excluded as pro- 
genitors for this event at the 90% confidence level. As a refer- 
ence point, a compact binary progenitor with masses 1 Mq < 
mi < 3 Mq and 1 Mq < m2 < 4 Mq with D < 3.5 Mpc is ex- 
cluded at 90% confidence; the same system with Z) < 8.8 Mpc 
is excluded at the 50% level. This result is averaged over dif- 
ferent theoretical waveform families; 20% of the simulated 
waveforms include spins sampled as described above. 

A number of systematic uncertainties enter into this analy- 
sis, but amplitude calibration error (^ 10%) and Monte-Carlo 
statistics have the largest effects. These uncertainties have 



been folded into our analysis in a manner similar to that de- 
scribed in (Abbott et al. 2005b, ). In particular, the amplitude 
calibration was taken into account by scaling the distance of 
the injection signal to be 1.28 x 10% larger; the Monte-Carlo 
error adds 1.28y^/7(1 -p)/n to p = p[0\h(t;m2^D)] where n 
is the total number of simulated signals in a particular mass- 
distance bin. 

We evaluate the hypothesis that the event occurred in M3 1 , 
as electromagnetic observations hint might be the case, given 
our observation. We adopt the measured distance of 0.77 Mpc 
to M31. We then simulated a large number of inspirals at 
distances 0.77 Mpc < D < 0.9 Mpc which allows us to ac- 
count for both uncertainty in distance to M31 (7%) (Freed- 
man et al. 2001) and the amplitude calibration uncertainty dis- 
cussed above. The simulations exclude any compact binary 
progenitor in our simulation space at the distance of M3 1 at 
the > 99% level. 



30 
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Fig. 3. — The probability as described in Eq. (3) is computed using in- 
jections made only into the 180 s segments immediately before and after the 
on-source time. The shaded regions represent 90%, 75%, 50%, and 25% ex- 
clusion regions, from darkest to lightest respectively. The distance to M31 is 
indicated by the horizontal line at Z) = 0.77 Mpc. Both amplitude calibration 
uncertainty and Monte-Carlo statistics are included in this result; apparent 
fluctuations as a function of mass are due to Monte Carlo uncertainty. 



4. SEARCH FOR A GRAVITATIONAL WAVE BURST 

To search for a gravitational wave burst associated with 
GRB 070201 we have used LIGO's current baseline method 
for near-real time searches for gravitational wave bursts as- 
sociated with GRB triggers (GCN 2007; IPN3 2007). A de- 
tailed description of the analysis method is presented else- 
where (Abbott et al. 2007a). 



4.1. Search Method 

The burst search method is based on cross-correlating a 
pair of pre-conditioned datastreams from two different grav- 
itational wave detectors. The pre-conditioning of the datas- 
treams consists of whitening, phase-calibration, and band- 
passing from 40 Hz to 2000 Hz. The cross-correlation is 
calculated for short time series of equal length taken from 
the datastreams of each detector. For discretely sampled 
time series and S2, each containing n elements, the cross- 
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correlation, cc, is defined as: 



^[^l(0-/il][^2(0-M2] 



\ 



(4) 



where fii and /i2 are the corresponding means of and 
S2. Possible values of this normalized cross-correlation 
range from -1 to +1, the minus sign corresponding to anti- 
correlation and the plus sign to correlation. 

The measurement of the cross-correlation statistic pro- 
ceeded as follows. Both 180 s, on-source time series of 
HI and H2 data were divided into time intervals (or cross- 
correlation windows) of length Tccw Previous analyses have 
shown that using two windows, Tccw = 25 ms and Tccw = 
100 ms, is sufficient to target short-duration signals lasting 
from ~ 1 ms to ~ 100 ms. The intervals were overlapped by 
half (i.e., Tccw/'^) to avoid missing a signal occurring near a 
boundary. The cross-correlation value, cc, was calculated for 
each H1-H2 interval pair and for both Tccw cross-correlation 
window lengths. The largest cc is the strength measure of the 
most significant correlated candidate value within the 180 s 
long on-source segment. To estimate the significance of this 
loudest event, we use off-source data to measure the cross- 
correlation distribution of the background noise. 

4.2. Background Estimation and Search Results 

Approximately 3 hours of data symmetrically distributed 
about the on-source segment were used to study the back- 
ground. These off-source data were collected from suffi- 
ciently close to the on-source time to accurately reflect the sta- 
tistical properties of the data within the on-source region. The 
detectors were collecting data continuously during the off- 
and on-source periods. The off- source data were divided into 
180 s long segments, corresponding to the length of the on- 
source segment. The off- source segments were treated identi- 
cally to the on-source segment. 

The distribution of largest cc values in the absence of a sig- 
nal was estimated for each cross-correlation window (Tccw = 
25 ms and Tccw =100 ms) by applying the method in Sec. 4. 1 
for all 180 s long off-source data segments. To increase the 
off-source distribution statistics, time shifts between the HI 
and H2 datastreams were also performed. The HI datastream 
was shifted by multiples of 180 s relative to H2. Then two 
180 s stretches from the two detectors were paired at each 
shift, making sure that two 180 s time stretches were paired 
only once. The distribution of cross-correlations obtained 
with time- shifted data is consistent with what is obtained from 
unshifted data. For both cross-correlation windows (Tccw), 
the resulting off- source loudest event cc distribution was used 
to estimate the probability that background noise alone (i.e., 
without a gravitational wave signal) would produce a cc value 
larger than the largest cross-correlation found in the on-source 
segment. 

Figure 4 shows the cumulative cross-correlation distribu- 
tion for the Tccw = 25 ms and Tccw = 100 ms cases. For the 
Tccw = 25 ms time- window, the largest cross-correlation found 
in the on-source data was cc = 0.36 (see arrow on Figure 4- 
a). The probability of obtaining a cross-correlation value this 
large from noise alone is 0.58. For the Tccw = 100 ms time- 
window, the largest cross-correlation found in the on-source 
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Fig. 4. — Cumulative distribution of measured cross-correlation values for 
the (a.) Tccw = 25 ms and (b.) Tccw = 100 ms cross-correlation windows. Both 
distributions with and without time shifts are shown, including the statistical 
errors. The arrows in both cases point to the largest cross-correlation found 
in the on-source segment. For the background distributions (1 -cumulative 
distribution) is plotted. 

data was cc = 0. 15 (see arrow on Figure 4-b), and the probabil- 
ity for this cross-correlation value is 0.96. These results are, 
therefore, consistent with noise. We conclude that no gravita- 
tional wave burst associated with GRB 070201 was detected 
by the search. 

4.3. Upper Limits on the Amplitude and Energy of 
Gravitational-Wave Transients Associated with 
GRB 070201 

Since the analysis of the previous section showed no evi- 
dence for a gravitational wave burst, we set upper limits on 
the amplitude and energy of gravitational waves incident on 
the detectors during GRB 070201. Denote the gravitational 
wave signal by h(t',hrss), where 



hrss=\ / {Mt)\^+\hy,(t)\^)dt 



(5) 



is the root- sum- squared amplitude of the gravitational wave 
signal. To determine an upper limit, one needs the probability 
of measuring cc given the presence of a signal with h^ss : 



p[cc\h(t;hrss)] 



(6) 



The search targets signals with duration < 100 ms. Within 
this class of signals, the sensitivity of the search has weak de- 
pendence on signal morphology; it depends primarily on the 
energy and the frequency content of the signal. Therefore, as 
long as the frequency and duration of the injected test wave- 
forms match the theoretical predictions, we can work with the 
waveform of our choice. A class of waveforms called sine- 
Gaussians have become the standard benchmark for burst 
searches and were used to construct the probability distribu- 
tion given in Eq. (6). The explicit formulae for h+(t) and /zx (0 
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are 



h+(t)= ho sin(27r/o(^-^o)) e^p 
hx(t)= ho cos(27Tfo(t-to)) exp 



-(27rfo(t-to)) 



21 



2Q^ 
-{2^foit-to)f 



222 



, (7) 
, (8) 



where fo is the central frequency, ho is the peak ampHtude 
of each polarization, is the peak time, and 2 is a dimen- 
sionless constant that represents roughly the number of cycles 
with which the waveform oscillates with more than half of the 
peak amplitude. Since the h+{t) and hx{t) waveforms have the 
same amplitude, these simulated gravitational wave bursts are 
circularly polarized. 

We provide results for the characteristic case of 2 = 8.9, as 
the dependence of the upper limits on Q is very weak. The 
measurement is carried out as follows. First, we choose a 
central frequency, /o, and an /Zj-ss value for the injected sig- 
nal. From these parameters, we calculate hit) using Eq. (1), 
Eq. (5), Eq. (7) and Eq. (8). We then add the calibrated h{t) 
to the on-source HI and H2 data, choosing a random start- 
ing time within the segments. We then measure the largest 
value of cross-correlation, cc, generally following the same 
method described in Sec. 4.1, except that in this case only a 
shorter interval around the injection is searched. Using the 
same /zrss values, we keep iterating the last two steps of the 
algorithm (randomizing a starting point and calculating the 
cc local maximum) until we have enough datapoints to deter- 
mine the conditional probability p{cc\hr^^). This probability, 
determined for different /zrss values and central frequencies, 
is then used to set a frequentist upper limit on /Zrss? given the 
largest cross-correlation found for the on-source segment in 
the search (see Sec. 4.1) (Abbott et al. 2007a). 

The resulting 90% /zrss upper limits are given in Table 1 
for circularly polarized sine-gaussians with different central 
frequencies and with Q = 8.9. The frequency dependence of 
the upper limits follows the shape of the detectorSs frequency 
dependent noise spectrum. 

The /zrss limits given in Table 1 include the calibration and 
statistical errors. These errors were propagated into the 90% 
^rss upper limits using the same procedure used in (Abbott 
et al. 2007a). The la errors considered were: (a.) calibration 
response phase error (10°); (b.) calibration response ampli- 
tude error (10%); and (c.) statistical error determined through 
Monte-Carlo simulations (2.1%). 

The upper limits on hrss implied by the burst search can 
be translated into conventional astrophysical units of energy 
emitted in gravitational waves. The gravitational wave en- 
ergy £^GW radiated by an isotropically emitting source that is 
dominated by emission at a frequency /o, is related to the /zrss 
received at distance D, much less than the Hubble distance, 
by (Riles 2004) 



771S0 , 



2 3 



2 ^2 ^2 



(9) 



Based on the sensitivity of this burst search as summarized 
in Table 1, we estimate that a gravitational wave burst with 
characteristic frequency in the most sensitive frequency re- 
gion of the LIGO detectors (/ ^ 150 Hz), if GRB 070201 
originated in M31 (at 770 kpc), must have emitted less than 
approximately 4.4x 1O~'^M0C^ (7.9x10^^ ergs) within any 
100 ms interval in the on-source window in gravitational 
waves. In terms of the SGR progenitor hypothesis, our ex- 
perimental upper limit on £'gw is several orders of magni- 



TABLE 1 

90% AMPLITUDE UPPER LIMITS AND CORRESPONDING 
CHARACTERISTIC ENERGIES FROM SINE-GAUSSIAN (SG) WAVEFORM 
SIMULATIONS ASSUMING 770 KPC AS SOURCE DISTANCE. THE hrss 
LIMITS GIVEN IN THE TABLE ALREADY INCLUDE THE CALIBRATION 
AND STATISTICAL ERRORS (ABBOTT ET AL. 2007a). 



SG central 
frequency (Hz) 



ms 



90% UL on 

hrss (Hz-1/2) 



100 
150 
250 
554 
1000 
1850 



25 
25 
25 
25 
25 
25 



2.15 X 10-2 
1.27 X 10-2 
1.34 X 10-2 
2.36 X 10-2 
4.12 X 10-2 
7.56 X 10-2 



100 
150 
250 
554 
1000 
1850 



100 
100 
100 
100 
100 
100 



1.97 X 10-2 
1.25 X 10-2 
1.41 X 10-2 
2.52 X 10-2 
4.51 X 10-2 
8.15 X 10-2 



Characteristic 

4;^ (Mqc2) 



Characteristic 



5.8 X 10-^ 
4.6 X 10-4 
1.4 X 10-3 
2.1 X 10-2 
2.1 X 10-1 
2.5 



(erg) 



4.9 X 10-4 

4.4 X 10-4 
1.6 X 10-3 

2.5 X 10-2 

2.6 X 10-1 
2.9 



1.0 X 10^1 
8.2 X 10^0 
2.5 X 10^1 
3.8 X 10^2 
3.8 X 10^3 
4.5 X 10^4 



8.8 X 10^0 

7.9 X 10^0 
2.9 X 10^1 
4.5 X 10^2 
4.7 X 10^3 
5.2 X 10^4 



o 



3 off-source 
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Fig. 5. — Differential (a.) and cumulative (b.) distributions of measured 
cross-correlation values for the Tccw = 100 ms cross-correlation window. Dis- 
tributions for both the 1.4M0-1.4M0 and 1.4M0-1OM0 binaries are shown 
along with the background distribution. The arrows points to the largest 
cross-correlation found in the on-source segment. On plot (b.) (1-cumulative 
distribution) of the off-source data is plotted. 

tude larger than the iC^^ergCD/TTOkpc)^ known to be emit- 
ted electromagnetically. And while present models for SGR 
bursts may differ substantially in their mechanism (de Freitas 
Pacheco 1998; loka 2001; Owen 2005; loka 2001; Horvath 
2005), they suggest that no more than 10"^^ erg is released in 
the form of gravitational waves. Therefore, the upper limit 
achievable with the present detectors does not exclude these 
models of SGRs at the M31 distance. 

We also estimate the sensitivity of the (100 ms) burst search 
to gravitational waves from a compact binary progenitor in 
M31 (see Figure 5). We choose as examples a IAMq-IAMq 
binary and a 1 AMq-IOMq binary. For each mass pair, we in- 
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ject approximately 1000 inspiral waveforms consistent with 
the distance of M31, with random isotropically distributed 
inclination and polarization, and with coalescence time uni- 
formly distributed through the on-source segment. Since, for 
these masses, the merger phase is expected to occur at fre- 
quencies well above that of maximum LIGO sensitivity, we 
inject only the inspiral portion. As for the sine-Gaussian sim- 
ulations, we determine the largest cross-correlation within a 
small time window around the coalescence time. We also 
account for possible systematic errors due to the calibration 
and the uncertainty in the distance to a possible source within 
M31, and statistical errors from the Monte Carlo procedure. 
We estimate with 90% confidence that a IAMq-IAMq bi- 
nary inspiral in M31 at the time of GRB 070201 would have 
a probability of at least 0.878 of producing a cross-correlation 
larger than the loudest on-source event. For IAMq-IOMq 
binaries this probability is at least 0.989. This gives us an in- 
dependent way to reject the hypothesis of a compact binary 
progenitor in M3 1 , while not relying on the detailed model of 
the inspiral signal. 

5. DISCUSSION 

We analyzed the data from the LIGO HI and H2 
gravitational-wave detectors, looking for signals associated 
with the electromagnetic event GRB 070201. No plausi- 
ble gravitational- wave signals were identified. Based on this 
search, a compact binary progenitor (neutron star + black hole 
or neutron star systems) of GRB 070201, with masses in the 
range 1 Mq < mi < 3 Mq and 1 Mq < m2 < 40 M©, located 
in M31 is excluded at the > 99% confidence. 

Our model-independent search did not find correlated sig- 
natures inconsistent with the noise within the HI and H2 data- 
streams that could be related to GRB 070201. Based on the 
sensitivity of our search and assuming isotropic gravitational- 
wave emission of the progenitor, an upper limit on the 
power emitted in gravitational waves by GRB 070201 was 
determined. A gravitational wave with characteristic fre- 
quency within the most sensitive range of the LIGO detectors 
(/ ^ 150 Hz) most probably emitted less than Egw <7.9x 
Iq50 gj^gg within any 100-ms-long time interval inside the on- 



source region if the source is in M3 1 . This limit on radiated 
power is comparable to the emitted power of some GRBs. 
However, it is significantly higher than the associated elec- 
tromagnetic emission of this particular GRB. Therefore the 
transient search only constrains the binary inspiral models for 
a short hard GRB in M3 1 and does not impose new limitations 
on magnetar-driven (SGR type) burst models. 

As gravitational- wave observations continue and the sensi- 
tivity of the instruments improves, we look forward to the as- 
trophysical insights that combined electromagnetic and grav- 
itational observing campaigns can bring. 
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